On the normalization of the HBT correlation function 
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We discuss the question of the normalization of the corre- 
lation function and its consistency with the often used form 
C(pi,P2) = 1 + |/(pi, P2)| 2 - We suggest an event mixing 
method which allows one to obtain absolutely normalized cor- 
relation functions from experimental data. 
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In the current note we address a question which being 
simple has caused confusing and erroneous statements 
in the literature. We discuss the question of the nor- 
malization of the correlation function used in the HBT 
analysis of multiparticle production. The normalization 
is, of course, a matter of convenience and author's favor. 
However, the question of the normalization becomes im- 
portant if one wants to compare the results from different 
papers or to compare the experimental data to the the- 
oretical predictions. In this note we show that in the 
literature there exist basically two definitions of the cor- 
relation function which differ by the normalization. Re- 
markably, starting from these two different definitions, 
authors arrive to the same result for the relation be- 
tween the correlation function and the particle emission 
(source) function. This means that some of the deriva- 
tions are wrong. 

The different definitions of the two-particle correlation 
function used in the theoretical papers on two-particle 
correlations in heavy ion collisions are the following. The 
first definition [0-0 is 
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(1) 
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where d 3 a/d 3 p and d e a/d 3 pi/d 3 p2 are the single-particle 
and the two-particle cross sections, respectively, and 
(7 is the total (inelastic) cross section, while d 3 n/d 3 p 
and d 6 n/d 3 pi/d 3 p2 are the particle and the pair yields 
per event, respectively. The relations between the 
cross sections and yields are d 3 n/d 3 p = l/a d 3 a/d 3 p, 
d 6 n/d 3 Pl /d 3 p 2 = l/a d 6 a/d 3 p 1 /d 3 p2- 

The second definition [&-pT| differs from the first one 
by the normalization factor: 
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C 2 (pi,p 2 ) 
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with (n) and (n(n — 1)} being the mean number of par- 
ticles and mean number of particle pairs per event, re- 
spectively 

The third definition is 
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where -P(pi) and P(pi, p 2 ) are called probabilities of ob- 
serving a single particle with momentum pi and a pair 
of particles with momenta pi and p 2 , respectively. Such 
a definition is not strict and may cause confusion. The 
word 'probability' can be understood in two ways: 1) 
P(pi)dpi is the probability of observing in a given event 
a particle with momentum between pi and pi + dpi ; in 
this case P(pi) = d 3 n/d 3 pi, or 2) P(pi) is the probabil- 
ity density for a given particle to have momentum pi . In 
the last case P(pi) = l/( n ) d 3 n/d 3 p±. Thus, depending 
on the interpretation of P(pi) and P(pi, p 2 ), the defini- 
tion (||) becomes equivalent to (|l|) or (||). In the rest of 
this note we consider the difference between (|l|) and <^). 

The factor, by which the definition (||) differs from 
||), is usually close to unity. For a Poissonian multi- 
plicity distribution it is exactly unity. In the presence 
of two-particle correlations, however, one cannot assume 
that the multiplicity distribution is Poissonian. Thus, 
in a general situation, the two definitions are different. 
Still, using these two different normalizations, different 
authors arrive to the same expression relating the corre- 
lation function to the source function. Brief review of the 
literature shows that the reason for this is that, in spite 
of declaring different normalizations at the beginning, ev- 
erybody (except, probably, Q where the normalization 
is kept through-out the paper) starts his/her calculation 
from the same expression for the correlation function in 
terms of the creation and annihilation operators (or anal- 
ogous one in terms of T-matrix) : 
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This equation, for chaotic source and Bose-Einstein 
statistics, leads (see for example |lj|) to the correlation 
function in the form: 
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C( Pl , P2 ) = l 
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= l + l/(pi,p 2 )| 2 (5) 



with the function /(pi,P2) depending on the shape of 
the boson source and asymptotically approaching zero 
for |p2 — pi | — *■ oo. For the Fermi-Dirac statistics there 
is a minus sign in front of the second term. 

The equivalence of the definitions (0) and (|^) can easily 
be shown using the expressions for one and two particle 
invariant distributions: 
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It also can be shown that the definition (||) and Eq. (||) 
are inconsistent. Indeed, combining them we obtain 
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Let us integrate this equation over pi and p 2 . Since the 
integral of d e n/d 3 pi/d 3 p2 is equal to the mean number 
of pairs per event (n(n — 1)), and the integral of d 3 n/d 3 p 
is the mean multiplicity (n), the integration of (Q) leads 
to 
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|/(pi, P2 )| 2 d 3 Pl d 3 p 2 = 0, 
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which is not possible unless /(pi, p 2 ) is identically equal 
to zero. Thus, the derivation starting from (||) and lead- 
ing to (jfy must contain an error. 

An expression analogous to (0) but obtained using 
the definition ([l]) would not contain the factor (n(n — 
l))/(n) 2 . An integration of this equation in absence of 
correlations gives (n(n— 1)) = (n) 2 . For identical bosons 
(without final state interactions) the correlation function 
is greater than unity, and thus the same integration yields 
(n(n — 1)) > (n) 2 which reflects the fact that bosons like 
to be produced copiously. For fermions the correspond- 
ing correlation function C(pi,p 2 ) < 1. This results in 
(n(n — 1)) < (n) 2 which reflects Pauli blocking. 

The whole presented problem is meaningless from the 
practical point of view as far as the experimentalists do 
not use Eq. (jl]) nor (||) to normalize the measured corre- 
lation functions. Instead, based on (|5|), they choose the 
normalization such that C(pi,p 2 ) — > 1 for |p 2 — Pi| — > 
oo. However, after we established that the definition 
(|l|) is correct in the sense that it matches (||), we may 
note that it does not leave any freedom of the normaliza- 
tion. For the often used event mixing technique it would 
mean the following. As usual, to generate mixed pairs 
one should take a pair of different events and combine all 
particles from one event with all particles from the other 
event. The normalization given by (|]J) will be automati- 
cally achieved if the number of event pairs used in event 
mixing is equal to the number of events used to calculate 
the actual pair distribution. 



Summarizing, we showed that the correlation function 
in the form (|^), often used in theoretical and experimen- 
tal papers, corresponds to the definition ([j]) and is not 
consistent with definition (||) . We suggest to use the def- 
inition (]l|) as the basis for the analysis of experimental 
data. 
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